Supersymmetric strictly isospectral FRW models for zero factor ordering 



H.C. Rosu and J. Socorro 
Instituto de Fisica de la Universidad de Guanajuato, Apdo Postal E-143, Leon, Gto, Mexico 



We apply the strictly isospectral technique of standard supersymmetric quantum mechanics to the Q — factor ordered 
Wheeler- De Witt equation for the closed Friedmann-Robertson- Walker (FRW) minisuperspace model. The resulting strictly 
isospectral one-parameter families of both FRW cosmological potentials and "wavefunctions of the universe" are exhibited with 
relevant plots. 
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Bianchi cosmologies |lj are cosmological models which 
are currently used to illustrate possible ways of finding 
exact solutions to the canonical constraints in quantum 
cosmology. In the last decade the number of exact solu- 
tions has been continuously increased by applying vari- 
ous formalisms || . The simplest of all cosmologies is the 
FRW model because of its full one-dimensional charac- 
ter. The goal of the present work is to apply the strictly 
isospectral scheme based on the general Riccati solution 
which is quite well known in nonrelativistic one- 
dimensional supersymmetric quantum mechanics, to the 
Q = factor-ordered Wheeler-DeWitt (WDW) equation 
corresponding to the closed FRW cosmological model. 
Thus, we will be able to display a one-parameter class 
of strictly isospectral cosmological FRW solutions and to 
present the wavefunctions of the universe for that case. 

FRW models belongs to the diagonal Bianchi type IX 
models in the class A models [Q. Therefore we recall 
briefly the canonical formulation in the Arnowitt-Deser- 
Misner (ADM) formalism of the latter ones. The metric 
has the form 



ds 2 = -dt 2 + c 2 ^) (e^y, 



(1) 



where a(t) is a scalar function and /3y(t) a 3x3 diagonal 
matrix, /3y = diag(x + \A3y,x — v3y, — 2x), lo 1 are one- 
forms which characterize each cosmological Bianchi type 
model, obeying du) 1 = §Cj k uP A w k , Cj k being the struc- 
ture constants of the corresponding invariance group. 
The ADM action has the form 



I 



P x dx + P y dy + P Q da - N7ij_)dt, 



(2) 



where the P's are canonically conjugated momenta, N is 
the lapse function and 



Ha 



-3a 



P 2 +P 2 +P 2 +e 4a V(x,y)) 



(3) 



c 4q V(x, y) = U(q M ) is the potential term of the cosmo- 
logical model under consideration. The WDW equation 
for these models is achieved by replacing P qf i by — id q ^ 



in Eq. (3), with q M = (a,x, y). The factor e _3a may be 
factor-ordered with P qf . in many ways. Hartle and Hawk- 
ing H have suggested an almost general factor ordering, 
which in this case would order e _3Q P 2 M as 

_ e -0-Q)o d aB -^d a = -e~ 3a d 2 a + Qe~ 3a d a , (4) 

where Q is any real constant. In the following, we shall 
work in the Q = factor ordering in which the WDW 
equation reads 

□ * - U(q M ) * = 0, (5) 
and by means of the Ansatz ^(q M ) = Ae ±$ one gets 

± Ad <f> + A[(V<I>) 2 — U] = 0, (6) 

where □ = G^^, (V) 2 = -(&)» + (f) 2 + (£)», 
and C" = diag(-l,l,l). 

For the FRW model there is only the a coordinate 
and the corresponding potential is U = — e Aa . Then the 
WDW equation can be written as a one-dimensional or- 
dinary differential equation 



da J \ da 2 



(7) 



and allows one to think of the Witten supersymmetric 
procedure (for a recent review see Q) since the change 
of function $ = ± In y produces a zero-energy one- 
dimensional Schrodinger equation, whereas the change 
^ = ±W leads to a Riccati equation for W, which can 
be identified with the superpotential in the Witten pro- 
cedure. The one-dimensional Schrodinger equation reads 



y 



e 4a y = 



(8) 



with the second derivative taken with respect to the 
a coordinate, and possessing as general solution a su- 
perposition of modified Bessel functions Zo(i^e 2a ) = 
C\Io{\e 2a ) + C2K (^e 2a ). Let us consider separately 
the two particular cases as follows: 
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(i) Zq = K ; This case corresponds asymptotically to 
an Ansatz of the type \& = Ae - * 1 for the "wavefunction 
of the universe", with $1 = \e 2a -Vol. Indeed, asymptot- 
ically (a » 1) § 



Ko(\e 2a ) 



— ~ e —a 



(9) 



(ii) Zq — Iq; Again asymptotically, it corresponds to 
an Ansatz of the type \I/ = Ae + ^ 2 for the 'wavefunction 
of the universe', with $2 = he 2a — a since for a 3> 1 @ 
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(10) 



We shall now apply the technique of strictly isospectral 
potentials to the FRW model. This technique is known 
since about a decade in one-dimensional supersymmet- 
ric quantum mechanics and usually requires nodcless, 
normalizable states of Schrodinger equation. However, 
Pappademos, Sukhatme, and Pagnamenta [^) showed 
that the strictly isospectral construction can be also 
performed on nonnormalizable states. Thus, using the 
nonnormalizable modified Bessel function Zq the strictly 
isospectral construction says that the isospectral poten- 
tials for the FRW cosmological model are given by 

S iso (a;X) =S-(a)-2[lii(J + X)]" 

4Z (a)Z' (a) 2Z 4 (a) 

= s {a) - Jo + x + (JoW' 

where S~ (a) — e 4a is the Schrodinger-FRW potential, A 
is a real parameter and 



Mn)= I Zl^-e^dy. 



(12) 



The common argument for the supersymmetric strictly 
isospectral method is as follows. One can factorize the 
FRW one-dimensional Schrodinger equation with the op- 
erators A = ^ + W(a) and A* = + W(a), where 
the superpotential function is given by W = — Z /Zq. 
Then, in the Witten scheme, the potential and superpo- 
tential enter an initial 'bosonic' Riccati equation S~ = 
W 2 — W . At the same time, one can build a 'fermionic' 
Riccati equation S + — W 2 + W , corresponding to a 
'fermionic' Schrodinger equation for which the opera- 
tors A and A + are applied in reversed order. Thus, the 
'fermionic' Schrodinger- FRW potential is found to be 

S + = e 4a — 2(-^-) . This potential does not have the 
Zq solution as a 'wavefunction of the universe'. How- 
ever, it is possible to reintroduce the Zq solution into the 
spectrum, by using the general superpotential solution 
of the 'fermionic' Riccati equation. The general Riccati 
solution reads 



The way to obtain Eq. (13) is well-known and will 
not be repeated here. From it one can easily get Eq. (11). 
The one-parameter family of 'wavefunctions of the uni- 
verse' differs from the initial one, being the A-depending 
quotient 



Z (i±e 2a ) 
Jo + X 



(14) 



da 



HM<*) + M 



(13) 



though in the limit A — > ±00 one recovers the initial solu- 
tion. In Eq. (14) we have included the factor A(A + 1), 
however not for normalization reasons as in supersym- 
metric quantum mechanics || (the wavefunctions are still 
nonnormalizable) but in order to discuss the so-called 
Pursey limit |hJ] A = 0, and the Abraham-Moses limit 
fllf A = — 1 (see figure captions). In order to get a better 
feeling of what is going on we have plotted the strictly 
isospectral FRW cosmological potentials and wavefunc- 
tions for Zq = Iq and Zq = Iq — Kq and two values of the 
A- parameter, see Figures 1, 2 and 3, confirming in this 
cosmological context some known facts of the supersym- 
metric strictly isospectral construction in nonrelativis- 
tic quantum mechanics What we actually did was 
to calculate the integral and all the other quantities in 
the interval (— c, a), where c is a sufficiently big positive 
number, while a has been chosen a not too big positive 
number. In other words, we have a particle (the FRW 
universe) in a special type of box, that we call a Bessel 
strictly isospectral box. For Iq 7 the integral Iq is a finite 
quantity even in the infinite limit, and the computer cal- 
culations show small well structures in the rising branch 
of the isospectral potentials and corresponding peaks in 
the wavefunctions. The Bessel box is single-walled at the 
a limit. For Jo — Kq, the integral is divergent in the infi- 
nite limit. In this case the Bessel strictly isospectral box 
is like having walls of different heights at the c limit and 
a limit. This peculiar box is a quite interesting object 
being the physical embodiment of the strictly isospectral 
family of FRW potentials as a family of Bessel solitonic 
potentials fl2|| possessing translationally moving struc- 
tures. It can also be considered as a kind of cavity with 
slowly moving walls and wells if A is varied even if the 
interval (— c, a) is kept fixed. 

In conclusion, we showed that the Hamiltonian con- 
straint of canonical quantum gravity for FRW minisu- 
perspace models may contain a one-parameter family of 
FRW cosmological potentials and "wavefunctions of the 
universe" . They are strictly isospectral (i.e. at zero en- 
ergy) WDW solutions of the FRW case including the 
common solution in the limit A — > ±00. Unfortunately, 
the strictly isospectral method can be applied only to ei- 
ther pure one-dimensional cases as the FRW model here 
(see also Q ) or to separable cases Q . 
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Fig. la 

The original potential and two members of the 
one-parameter strictly isospectral family of FRW potentials 
in the case Zq — Io for A = 25 (small well) and 100 (bigger 

well) 
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Fig. lb 

The corresponding wavefunctions. 



Fig. 2a 

The original potential and two members of the 
one-parameter strictly isospectral family of FRW potentials 
for To — Ko, and the same values of the A-parameter. On the 
arbitrary vertical scale the two peaks close at a value of 
about 2,000. 
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Fig. 2b 

All three potentials on the same plot, showing at the left the 
peak corresponding to A = 100 and slightly at the right that 
for A = 25. We have found that the peaks are moving to the 
left if A is positive and increasing (i.e., beyond Pursey limit), 
while they are moving toward the origin if A is less than the 
Abraham-Moses limit and decreasing. 



Fig. 2c 

A more detailed plot of the region close to the origin for the 

three potentials showing the real trend of the plateau 
between the small well structures at the right and the peaks 
at the left. 
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Fig. 3 

The strictly isospectral 'wavefunction of the universe', Z gm , 
for Zq = 7o — Ko and the A parameters as in Figs. 1 and 2 
(the small peak corresponds to A = 25, whereas the big one 
is for A = 100). We notice also that the positive tails (at the 
right of the peaks) are practically zero, while the negative 
tails have a small positive slope. That implies that the 
wavefunctions are not normalizable. The segment very close 
to the vertical axis is a portion of the original 
'wavefunction', close to the value Iq(1/2) — Kq{1/2). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9606030v4 



This figure "figl-2.png" is available in "png" format from: 
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